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SUMMARY 

An analysis of low-speed boundary— layer flow over airfoils 
at moderate angles of attack and at Eeynolds numbers of several 
million is presented. Methods are developed for estimating the 
growth of the boundary layer In the regions of laminar separation 
and transition. Calculation of the growth of the turbulent 
boundary layer and the chordwise location of the turbulent 
separation point are considered. 

The conceptp, which are found to be basic for unders tending 
low-speed boundary— layer flow, are also applied to data for 
transonic speeds. A close similarity Is found to exist between 
low-speed and transonic boundary— layer flow. An approximate 
procedure for calculating boundary— layer growth through shock 
waves is presented. The local effects of both laminar and 
turbulent boundary— layer separation are considered. 


INTRODUCTION 

Experimental measurements at transonic speeds have indicated 
that changes In the airfoil surface condition or of Eeynolds 
number cause marked changes in the chordwise static— pressure 
distribution. In particular > the flow in the vicinity of the 
shock wave is entirely different for laminar than for turbulent 
boundary— layer flow. Before an adequate analysis of transonic 
flow past airfoil sections can be made, it is necessary to obtain 
same understanding of the characteristics of boundary layers at 
transonic speeds. 

There Is little analysis available for boundary layers at 
transonic speeds; however, a vast literature dealing with low— speed 
boundary— layer theory exists. The question naturally arises as to 
whether some, of these numerous methods of analysis can be extended 
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to apply to transonic flow. The various houndary-layer flow regimes 
at low speeds are considered for the case of an airfoil at a 
moderate angle of attack. A number of basic characteristics of 
boundary— layer flow appear in this analysis. It is the purpose of 
this report to examine some transonic boundary— layer data to 
determine whether fundamental similarity exists between boundary- 
layer flow at low and transonic speedy. 

An airfoil at a moderate angle of attack was selected for 
the low— speed case in order to avoid the complex problem of 
predicting the chordwise location of transition from laminar to 
turbulent flow when this occurs ahead of the theoretical position 
of laminar separation. In the case considered, transition from 
laminar to turbulent flow follows separation of the laminar 
boundary layer and subsequent reattachment of the flow to the 
airfoil surface. This case is fruitful because regions of both 
laminar and turbulent separation occur. 

There are several semiempirical methods for calculating the 
development of the turbulent boundary layer and the position of 
turbulent separation. Data obtained at the Bureau of Standards 
and recently presented by Dryden (reference 1) assist-in evaluating 
the methods for calculating turbulent separation. These data show 
that the skin— friction coefficient of turbulent boundary layers is 
highly dependent upon the velocity profile, that is, the shape of 
the velocity distribution through the boundary layer. This 
dependence was not established at - the time von Doenhoff and Tetervin 
(reference' 2) developed their semiempirical equation for calculating 
the chordwise position of turbulent-separation. It appears that the 
complexity of their equation is at least partly due to attempting 
to fit an equation to the experimental data while neglecting the 
dependence of skin— friction coefficient on the velocity profile. 

The much simpler procedure developed by Gruschwitz (reference 3) 
is in substantial agreement with the data of reference 1 and Is a 
fair approximation to the data of reference 2. Therefore, the 
semiempirical equation of Gruschwitz Is used In the present report’ 
for calculating the turbulent separation point. 

A significant experimental study of boundary— layer flow at 
transonic speeds has been made by Ackeret, Feldmann, and Bott 
(reference h. ) These data are considered in the present report 
In terms of the concepts developed for low-speed flows. 
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airfoil chord length 
surface bump height 
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ratio of boundary— layer displacement to momentum thickness 

(s*/e) 

length of region of laminar separation 
ratio of local velocity to local velocity of sound 
ratio of local velocity to critical velocity of sound 
static pressure 
total head 
dynamic pressure 

Reynolds number based on airfoil chord 
Reynolds number based on boundary— layer momentum thickness 



pressure coefficient 

local velocity in boundary layer 
local velocity outside boundary layer 
free— stream, velocity 
maximum perturbation velocity 
chordwise distance 
distance normal to surface 

velocity gradient 

total boundary— layer thickness 
boundary— layer displacement thickness 

[•‘-.i’ 8 C 1 -^) ■*] 

density outside boundary layer 
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P u density in boundary layer at point where velocity equals u 

6 boundary— layer momentum thickness 

{■■£ 

T) G-ruschwitz boundary— layer shape parameter (equation (6) ) 

t 0 surface shearing stress 

t q / 2 q skin^-friction coefficient 

V kinematic viscosity 


AWALYSIS OS’ LOW— SEEED BOUWDAEY-LA.YER FLOW 

Consider a smooth— surf ace d airfoil at a moderate angle of 
attack in a low— turbulence low— speed flow. The various regimes of 
boundary— layer flow over the airfoil can be readily differentiated 
on the basis of the local pressure distribution. To avoid the 
confusion of plus and minus signs, the pressure coefficient S, which 
is defined as the ratio of the local to the free— stream velocity 
squared, will be used. The maximum pressure coefficient occurs on 
the upper surface near the airfoil leading edge. (See fig. 1.) 

Over the next few percent of the airfoil chord, there is an abrupt 
decrease in S. This decrease occurs in two parts which are divided 
by a short region of constant pressure. The corresponding behavior 
of the boundary layer is as follows : From the peak pressure point 

to the constant— pressure region the flow is laminar. In the region 
of constant pressure, which will be shown to be about - ' 1 percent of 
the chord in length for Eeynolds numbers of about one million, the 
boundary layer remains laminar but is separated from the surface by 
a " dead-air" bubble. The flow then turns turbulent and spreads 
back to the surface with the accompanying second abrupt decrease 
in S. Aft of the sharp forward pressure peak there are relatively 
moderate chordwise pressure gradients (fig. l) and the boundary 
layer is turbulent. Over the rear portion of the airfoil there 
is a region of relatively constant pressure. In such a region, S 
is somewhat greater than one, and the local boundary— layer flow is 
referred to as turbulent separation. 

The pressure distribution which has just - been described is 
the actual viscous— flow pressure distribution which approximates 
the potential— theory pressure distribution exoept in the vicinity 
of both the region of laminar and turbulent separation. 


The various regimes of flow will now be treated in detail 
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Laminar Flow 

The extent of. the laminar flow on the airfoil far the case 
"being considered is small but it is of importance for determining 
the initial conditions in the turbulent boundary layer. In refer- 
ence 5, von Barman and Millikan have studied laminar boundary— layer 
separation for the case of pressure distributions having a sharp 
peak. It is shown In this reference that the conditions for 
laminar separation are independent of both Eeynolds n umb er and the 
actual m ag n itude of the peak pressure coefficient. When the maxlTmiTn 
pressure occurs near the leading edge, as in the present case, the 
laminar bou nd a r y layer separates at the chordwise station at which 
the local value of S is about 0.8l times the peak pressure coeffi- 
cient. ' Furthermore, as can be determined from the results of 
reference 5, the local boundary— layer momentum thickness at this 
separation point is simply 


e 0.4 

c 


( 1 ) 


where c is the chord length, R Q Is the airfoil Eeynolds 
number, and p is defined as 

d «/~S _ d(u/tT 0 ) 

P d(x/o) d(x/c) 

For the low— speed type of airfoil pressure distribution being 
analyzed in the present report, the value of p is nearly constant 
from the pressure peak to the laminar separation point so that 
equation (1) is applicable. From reference 5, it is also found 
that at the laminar separation point the ratio of boundary— layer 
displacement to momentum thickness has a value of 3.8, 

A region of separated laminar flow Is characterized by 
constant local pres sire . Thus, as there are no pressure or shear 
forces acting on the boundary layer, there cannot be any appreciable 
changes In the boundary— layer momentum defect; that is, the 
boundary— layer momentum thickness must remain virtually constant 
over the constant— pressure region. At the same time there must 
be a rapid local increase in the boundary— layer displacement 
thickness since It is necessary to have a distortion of the 
streamlines to decrease the local pressure gradient to zero. In 
reference 6, laminar separation was measured on an RACA 66,2— 216 
airfoil section at an angle of attack of 10.1° and for a Reynolds 
number range of 0.9 to 2.6 million. An analysis of these data 
Indicated that the momentum thickness does remain virtually constant 
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in the region of constant pressure vhile there is a marked, increase 
in the "boundary— layer displacement thickness. (See fig. 2 .) 


There appears to "be the following factor determining the 
magnitude of this Increase in the displacement thickness. Through- 
out the range of Eeynolds numbers for which data are presented, the 
length of the region of constant pressure can "be characterized by 
a value of Rl equal to approximately 25,000, where R^ is the 
length of the region of constant pressure times the local velocity 
and divided by the kinematic viscosity. Two other properties of 
the laminar separated region are known; The local pressure, hence 
velocity, are constant; whereas in the absence of separation there 
would be a large local velocity gradient p. The shape and 
magnitude of the dead-air bubble are such as to bring about this 
change in the local velocity gradient. From thin airfoil theory, 
it is known that a local velooity perturbation is proportional to 
the height of a local bump 

( 2 ) 


h _ ZfiJ 
L~ Kr U 


where 


h height of bump 

NJ velocity perturbation due to bump 

U local stream velocity 

L length of region of constant velocity characterized by 

R l - 25,000 

The following data obtained from reference 6 indicate that the 
value of K is about one : 


Reynolds 

number 

U 

h 

L 

0.9 X 106 

1.5 x 10 s 

2.2 X 10 6 

0.10 

.10 

.08 

0.08 

.08 

.07 


The height of the dead-air bubble is thus directly determinable 
from the potential— theory pressure distribution and equation (2) . 
That this height is equal to the increment of the displacement 
thickness over the constant— pressure region follows directly from 
the definition of the displacement thickness, since the air velocity 
inside the bubble is virtually zero. 
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The boundary— layer momentum and displacement thickness at the 
end of the region of 3 wmi t-iat" separation are thus calculable . The 
looal momentum thickness is obtained from equation (1), and the 
displacement thickness is 3-8 times this value pluB the maximum 
displacement of the separated flow from the surface . The latter 
quantity is obtained from equation (2) . 


Transition Region 

At the termination of the short constant-^ pres sure region, the 
separated laminar flow becomes turbulent and spreads back toward 
the airfoil surface . Over the region in which the separated 
turbulent flow spreads toward the airfoil surface to establish 
itself as a turbulent boundary layer, there is a large pressure 
gradient. An analysis of the experimental boundary— layer profiles 
in reference 6 indicated that over this region there was only a 
small variation in the boundary— layer displacement thickness. This 
variation consisted of an initial decrease followed by an increase, 
so that the displacement thickness at the end of this region of 
large pressure gradient was nearly the same as at the beginning. 

At the same time, the boundary— layer momentum thickness increased 
rapidly. The classical equation for analyzing the growth of 
momentum thickness is the boundary— layer momentum integral equa- 
tion 

ds (se+6*) as _ jo . . 

dx 2S dx 2q W 

where 


e 

momentum thickness 

X 

length along airfoil surface 

8* 

displacement thickness 

S 

pressure coefficient 

1 

local dynamic pressure 

T o/2q. 

skin— friction coefficient 


The growth of the boundary— layer momentum thickness in the transi- 
tion region can be calculated by means of equation ( 3 ) . This 
calculation is simple because when the flow is detached from, the 
surface the skin— friction coefficient is zero; moreover the experi- 
mental data suggest the use of a constant value of the boundary— layer 
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displacement thickness, namely, its value at the end of the region 
of la m inar separation. Hence equation (3) reduces to 



where subscripts i and s indicate, respectively, conditions at the 
beginning of and at some later point in the transition region. 

In figure 2, a comparison is presented between the theoretical 
and experimental momentum, thickness distribution in the region of 
laminar separation and transition for the KACA 66,2— 216, a = 0.6 
airfoil section at an angle of attack of 10.1° and a Reynolds number 
of 900,000. According to the theory, the momentum thickness should 
remain constant over the region of separation at a value given by 
equation (1) . In the transition region the momentum thickness growth 
Vas calculated by means of equation (If), using the experimental 
pressure distribution obtained from reference 6 and assuming a 
constant value of displacement thickness of 0.0009 chord lengths. 

The value of momentum thickness at the point of laminar 
separation estimated, using equation (1), agrees well with the 
experimental value. However, there is a growth of momentum thickness 
in the region of constant pressure not predicted by the theory. 
Consequently, the theoretical values of momentum thickness in the 
following transition region are lower than the experimental because 
they depend upon the theoretical value at the end of the constant- 
pressure region which is somewhat lower than the experimental. 


Turbulent Flow 

A number of methods are available for calculating the growth 
of the turbulent boundary layer. These methods consist of different 
procedures for solving the boundary— layer momentum equation. The 
factors which appear in this equation are the chordwise pressure 
distribution, the local ratio of displacement to momentum, thickness, 
and the looal skin-friction coefficient. At points some distance 
ahead of the turbulent separation there is uncertainty in the 
theoretical estimation of each of these faotors. The presence of 
a region of separated flow appreciably alters the pressure distri- 
bution, ahead of as well as behind the separation point, from the 
unseparated or potential— flow distribution. The ratio of displace- 
ment to momentum thickness varies from about 1.4 to 2.5 over the 
region of turbulent flow. The data recently presented by Dry den 
show that the variation of this ratio results in a marked variation 
in the skin-friction coefficient. 
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In spite of the previously mentioned difficulties, relatively 
simple approximate solutions of the momentum equation can he used 
to obtain quite accurate values for the boundary— layer growth up 
to the vicinity of the turbulent separation point. In cases where 
turbulent separation occurs, errors in evaluating the skin— friction 
coefficient are of secondary importance because the effects of the 
local pressure forces on the boundary— layer growth are larger than 
the surface shear forces. The ratio of displacement to momentum, 
thickness 6*/0 appears in the momentum equation only in the group 
(6 */6 + 2) which, although ranging from 3.^ to k. 5 , only varies 
between 3 A and 3*8 for the major portion of the region of 
unseparated turbulent .flow. It thus appears that the greatest 
uncertainty in the theoretical calculation of boundary— layer growth 
arises from changes in the pressure distribution due to the presence 
of a region of turbulent separated flow. The magnitude of these 
pressure changes depends upon the variation of the boundary— layer 
displacement thickness in the wake as well as over the rear portion 
of the airfoil. This is a complex problem which has not as yet 
been adequately solved. Of course, if the boundary— layer calcula- 
tions are based on an experimental pressure distribution, this 
complexity is circumvented. 


In view of the uncertainty in evaluating the basic parameters, 
it follows that elaborate step— by— step calculations of the boundary- 
layer growth involve unwarranted complexity. For fully developed, 
turbulent, boundary— layer flow, a reasonable value for the skin- 
friction coefficient can be obtained from Falkner's equation for 
turbulent flow over flat plates (reference 7)/ namely 


2<1 153 Be 1/8 

where Eg is the local Beynolds number in which the momentum 
thickness is used as the characteristic length. The chordwise 
variation of boundary— layer momentum thickness can be obtained 
most directly by substituting this equation into equation (3) and 
integrating under the assumption of a constant average value of H 
of 1.6 which leads to 



O.OO76 



2 

S 2 d (x/c) + 



( 5 ) 


A knowledge of the chordwise growth of boundary— layer momentum 
thickness gives no indication of separation. Boundary— layer 
measurements in the vicinity of the turbulent separation point are 
in general agreement that the boundary— layer shape parameter 
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(H = 6 */S) can be used as an Indicator of the location of the 
turbulent separation point. The value of H at separation is 
between 1.8 and 2.7, usually being about 2.5. Therefore, a 
relationship is needed to permit the calculation of the boundary- 
layer shape parameter. In reference 2, the semiempirlcal equation 
derived for calculating the variation of H Is 


e 


dH 

dx 


4.680(H-e.S7s) 

e 


0 dS 2q 


2.035 (H - 1.286) 


where — is the equivalent flat-plate skin— friction coefficient 
2q. 

based on the local value of Bg. It appears that this equation is 
excessively complex as a result of the fact that, in fitting an 
equation to the experimental data, it was assumed that the value 
of the skin-friction coefficient was independent of H. The data 
presented in reference 1 indicate that, for values of H from about 
1.5 to 2. 2, the skin-friction coefficient is approximately 

proportional to ( T 0 /2q) / (E-1.29) • Thus, when the dependence of 


skin-friction coefficient on H is taken into account the von Doenhoff— 
Tetervin equation may reduce to the much simpler form 


F(H)dH = — ^ — A — 

S t o 0 

However, there are not sufficient data available to permit an 
accurate determination of the variation of skirt-friction coefficient 
with H as well as with Reynolds number. In the present report 
the Grusohwitz equation will be used. 

Gruschwitz assumed that the turbulent boundary— layer velocity 
profiles formed a one— parameter family. The shape parameter he 
used was defined as 

n = i-(g) (6) 

\ TJs y=0 

This parameter is related to 0 , x, and S by the empirical equation 


0 d(Srj) 


0.00461 - 0.00894 TJ 


S dx 
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The data of references 1 and. 2 lend, considerable weight to the basic 
assumption of Gruschwitz. In both of these references, the velocity 
distribution through the boundary layer is shown to depend only 
upon H. The variation of the boundary— layer velocity profile with 
H is studied by plotting curves of u/tJ versus H for various 
constant values of y/©. The data in these two reports are in good 
agreement for values of jfo above one. 

An analysis of the data of references 1 and 2 revealed that 
for each airfoil the quantity St} varied less than 20 percent 
throughout the region of turbulent flow. This variation appeared 
to be somewhat dependent upon Reynolds number and hence these data 
were not in close agreement with the equation presented by 
Gruschwitz. However, for Reynolds numbers of several million the 
average variation for each configuration investigated was well 
represented by 


9 d(STj) . v 

= 0.005 - 0.009 n (7) 

b CLX 

which is essentially the G-ruschwitz equation. This equation was 
investigated Only for regions of decreasing pressure coefficient S. 

It can be used to calculate the turbulent boundary— layer separation 
point for airfoils at moderate angles of attack. The value of rj 
for which fully turbulent flow starts is approximately O.58 (H = l.h), 
and separation occurs when tj attains a value of. about 0.93 (H= 2 . 5 ). 
In- performing this computation, it is merely necessary to know the 
chordwise pressure distribution. W)ien the region of separated 
turbulent flow is more than 10 percent of the chord In length, there 
is a substantial difference between the actual and the potential- 
theory pressure distributions. Consequently, if the calculations 
are based on the potential— theory pressure distribution, considerable 
error can exist in estimating the chordwise position of the separa- 
tion point. However, since, as was noted above, there is only a 
moderate chordwise variation in the quantity Stj, it Is possible 
to obtain a close estimate of the value of the pressure coefficient 
S over the separated region. In fact, a useful rough approximation 
is 


(an) 


*1=0.58 



*1=0.93 


or 


%= 1.4 2 ! 2 ~ 

S H=2.5 U2 E=2.5 
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which emphasizes that turbulent separation. Just like laminar 
separation, is primarily a result of decreases in the pressure 
coefficient S. 

The ohardwise variation of tj (or H) is of interest for 
evaluating cha ng es in the velocity profile shape as well as for 
predicting separation. An interesting approximation far these 
velocity profiles can be obtained as follows: 

Assume that the velocity distribution through the boundary 
layer is related by a power law to the distanoe from the surfaoe, 
namely. 



From, the 


definitions of the various quantities it 

8 - V _ N 

8 = (1+N) (1+2N) } 8 = Jhl ’ 


follows that 
1 + 2N 


and hence 


r r iiif 

u 7 g L 

U'l 0 L (l+N) (1+2N) J J 


(8) 


This simple expression is in surprisingly oloBe agreement with the 
experimental data. The G-ruschwitz shape parameter is then 


n 



= i - 


ar 

(l+N) (1+2N) 



(9) 


The variation of tj with H given by this expression is shown in 
figure 3, which agrees satisfactorily with the curve obtained 
empirically by Grusohwitz. A further indication of the value of 
using ij as a parameter is that the variation of skin-friction 
coefficient with ij is linear for the data presented by Dryden. 
These data indicate a relation between the skin-friction coeffi- 
cient and the boundary— layer shape factor such as 

£ - £ < 2 - 8 - *•> 

T 

where — ^ is the flat-plate skin-friotion coefficient. 

2q 


Hie turbulent separation on an NACA 66,2-216, a = 0.6 airfoil 
section at 10.1° angle of attack and a Reynolds number of 2.6 million 
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is a n a lyzed in reference 2. The experimental pressure distribution, 
for this airfoil is shown in figure 1. The growth of the momentum 
thickness in the turbulent region is computed by means of e ques- 
tion ( 5 ) and the experimental pressure distribution. The values ~ 
obtained in this computation are substituted in equation ( 7 ) and the 
chordwise variation of tj and hence H evaluated. These values 
of H and also the calculated values of momentum thickness are 
compared with the experimental data in figure h. Although the 
computed values of H are smaller than those measured experimentally, 
the curves are similar in shape, and the fact that turbulent separa- 
tion is imminent at the 0.7 chord station is predicted by the theory. 


ANALOGIES BETWEEN BOUNDARY— LAYER FLOW AT 
TRANSONIC AND AT LOW SPEEDS 

A number of the concepts which were developed in the preceding 
section for low— speed flows will be shown to apply in compressible 
flow. Recently Aokeret, Feldmann, and Rott (reference 4) published 
a thorough experimental investigation of boundary— layer flow in the 
vicinity of compression shock waves. These data will now be examined 
in the light of the preceding low— speed analysis. At transonic speeds 
the flow over an airfoil may be divided into three regions : the 

forward subsonic region, followed by a supersonic region terminated 
by a compression shock which returns the local flow to subsonic 
speeds. It Is the character of the return of the flow from super- 
sonic to subsonic velooity which varies with the state of flow in 
the boundary layer. 


laminar Plow 

Consider first the case of laminar boundary— layer separation 
In the vicinity of a compression shock wave. When the boundary 
layer Immediately ahead of the shock wave Is laminar and the local 
Mach number is about 1.2 or greater, the compression shock wave 
Is shaped like the Greek letter X. The front leg of the X is 
an oblique shock wave arising from the deflection of the boundary 
layer from the surface, that Is, 1 ami nar separation. The rear leg 
is the main compression wave through which the flow goes from 
supersonic to subsonic velocity. Between these two branches of 
the shock wave the surface static pressure remains constant, a 
characteristic property of regions of laminar separation. Ackeret, 
Feld.mann, and Rott present boundary— layer measurements for this 
type of flow. No measurenents of the local velocity within the 
bubble of separation are presented. They assume appreciable 
negative velocity to exist in this bubble; whereas these velocities 
should be considered negligible. The boundary-^ lays r— momentum. and 
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displacement-thickness distributions were recalculated for zero 
velocity in the dead— air bubble. In performing these recalculations, 
the Acker et procedure was used for obtaining the approximate density 
variation through the boundary layer from the measured Mach number 
j distribution. This procedure is based on the assumption that energy 
is constant through the boundary layer. These revised experimental 
values are shown in figure 5. It is seen that over - the constant- 
/ pressure region the displacement thickness grows approximately 
linearly (a phenomenon to be discussed later), while the momentum 
thickness remains virtually constant. Over the transition region 
of rapid pressure recovery, at the base of the main shock wave, 
displacement thickness remains constant but momentum thickness 
increases. The boundary— layer behavior In these regions is similar 
to that previously noted for low— speed laminar separation and 
transition. Moreover, at Mach numbers of about 1.2 the length of 
the constant-pressure region can be characterized by a Reynolds 
1 number run of 100 , OCX), a value four times that observed at low 
J speeds. This particular numerical value is characteristic of 
every case of transonic laminar separation presented in reference k. 
The difference between the length of the laminar separated region 
at low and high speeds may be due to the increase in the stability 
of laminar flow with Mach number (reference 8) . 

Another analogy between transonic and low— speed laminar 

separation Is the local flow deceleration wjbich is associated with 

the onset of laminar separation. Von Karman and Millikan have 

shown that at low speeds laminar separation occurs at the point 

where U 2 . . Aj 2 . attains a certain numerical value, 

separation' maximum 

This value depends only upon the chordwise position of the maximum 
local velocity and, for cases in which this position is well 
removed from the leading edge, the numerical value of the ratio is 
approximately 0 . 88 . This analysis could be expected to apply at 
transonic speeds because U 2 /2 is the kinetic energy per unit mass 
of fluid in compressible as well as incompressible flow. The 
numerical value of this ratio for the transonic pressure distribu- 
tion shown in figure 5 is 0.92 which Is in satisfactory agreement 
with the low— speed analysis. Moreover, this ratio corresponds to 
a flow deflection of 1.5° when supersonic oblique— shock-wave 
theory is applied. In figure 5, this value is presented an the 
theoretical curve for the growth of the boundary— layer displacement 
thickness over the region of laminar separation. The theoretical 
variation of the displacement thickness is seen to be only about 
one-half of that found experimentally. This difference is primarily 
due to the curvature of the surface. The flow is deflected 1 . 5 C 
with reference to a line tangent to the surface a-t the point of 
deflection. Consequently, In addition to the increase in displace- 
ment thickness caused by the flow deflection, there is an increase 
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measured "by the perpendicular distance from the surface to the line 
tangent to the surface at the deflection point. It was not possible 
to determine this distance from the data presented in reference ij-. 

A more exact calculation would consider this effect of surface 
curvature and then assume that 5* is constant over the following 
transition region. The momentum thickness is theoretically constant 
over the region" of laminar separation while the experimental measure- 
ments show a small increase. In the region of transition and 
reattachment , which occurs at the base of the main shock wave, the 
transonic data indicate that the displacement thickness is virtually 
constant. This result Is the same. as that noted for low— speed flow. 

In order to calculate the growth of the boundary— layer momentum 
thickness in the transition region. It is necessary to employ the 
momentum integral equation. For a compressible fluid, this equa- 
tion Is 


dS 9 dp 26+6* cU 2 t q 

dx p dx 2U S dx 2q 


( 10 ) 


where p and U are, respectively, the local density and velocity 
Immediately outside the boundary layer. In the transition region, 
the flow is detached from the surface so the skin-friction coeffi- 
cient Is zero. As a first approximation to the density variation 
outside the boundary layer adiabatic variation is assumed. Then 
the solution of equation (10) is 

^ /K _ Ms* 

'8 1 *+K8 s S Mi* 


where subscripts 1 and a are, respectively, conditions at the 
beginning and end of an interval in which K is the average value 
of 


12 — 7 M * 2 
Z " S-M*2 ~ 

The momentum-thickness distribution calculated, by means of this 
approximate solution, over the transition region is also presented 
in figure 5* 

A further effect of laminar separation, which Is indicated by 
Ackeret f s data, is a rearward movement of the position of the main 
shock wave from Its position when the houndary layer is turbulent. 
Consider the subsonic chordwise velocity distribution behind the 
shock wave in the latter case as a reference distribution; it is 
apparent that the presence of the effective bump, which is the bubble 
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of laminar separation, raises the local velocities above the 
reference distribution. Thus, the ohordwise location of the point 
at which sonic velocity occurs moves aft, and hence the shock wave 
moves aft when fully turbulent boundary— layer flow is replaced by 
laminar separation. It should be noted that between these two types 
of shocIMrave b oundary— layer interaction there is another possible 
type. When laminar separation occurs some distance from the leading 
edge of an airfoil, the laminar boundary layer at the separation 
point is rather thick and easily destabilized. Thus, the very process 
of laminar separation could precipitate transition ahead of the main 
shock wave. This would result in a marked change in the appearance 
of the main shock wave from that of the lambda-type wave. 

The question arises as to whether the boundary— layer 
momentum integral equation can be used to calculate the boundary- 
layer growth through a shook wave. The primary factor which could 
invalidate the use of this equation would be the occurrence of 
static pressure variation through the boundary layer. In refer- 
ence 4, static pressure surveys through the boundary layer are 
presented. The experimental measurements show that the static 
pressure is virtually constant through the boundary layer even 
when there is a large static pressure gradient Immediately outside 
the boundary layer. Equation (10) is therefore applicable. However, 
before attempting a solution of this equation for turbulent flows, 
it is necessary to investigate the characteristics of turbulent 
boundary layers at high speeds. 


Turbulent Flow 

The transonic, turbule nt-^b oundary— layer velocity profiles 
presented by Ackeret can be analyzed by the methods used for low- 
speed turbulent velocity profiles. The variable considered was 
the ratio of the local velocity at various points inside the 
boundary layer to the local velocity immediately outside the 
boundary layer. The points for which the local velocity ratio 
was considered were at a distance y from the airfoil surface 
equal to the local momentum thickness 0 , and also at even integral 
multiples of this distance y/0 from the surface. The latter points 
are plotted against the former in figure 6 . The curves drawn in 
this figure were obtained for lew— speed turbulent boundary layers. 

The data for a value of y/0 of, 4 indicate that the boundary- 
layer velocity profiles at transonic speeds may differ somewhat 
from those at low speed over the outer portion of the profiles. 
However, at the inner portion of the boundary layer, the values for 
transonic speeds are in agreement with the curves obtained from 
low— speed data. It therefore follows that essentially the same 
boundary— layer velocity profiles occur at low and transonic speeds, 
and throughout this Mach number range the same Gruachwitz shape 
parameter tj applies. 
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The relation between 'and 5*/e (or H) is different for 
c empress it le and incompressible flows because of the variation of 
density through the boundary layer. Ache ret assumes that the 
total energy and static pressure are constant through the boundary 
layer. From this it follows that the density variation through 
the boundary layer is a function of the velocity— profile shape, 
the free— stream Mach n umb er, and the local stream Mach number M. 
Values of H have been calculated for free— stream Mach numbers 
between 0.7 ami 1.0, and various local Mach numbers. As long as 
the local Mach number is greater than the free— stream Mach number 
but less than 1.4, the effect of compressibility on the value of 
H is well approximated by 


H = B*t-o(l + | M 2 ^) (11) 

It is now possible to consider solutions of equation (10) which 
can be used to calculate the turbulent-b oundary-lnyer growth through 
a shock wave. The pressure and Mach number variation in this region 
are large, anti, both of these factors have appreciable effects on H. 
The pressure variation causes an increase in H, the magnitude of 
which is reduced as a result of the decreasing local Mach numbers. 
This suggests the simple approximation of again using a constant 
average value of I in solving the momentum integral equation. This 
is such a gross approximation that it would be foolish to attempt 
to include the relatively small effect of the surface shear forces. 
Assuming a constant average value of 6 */0 and neglecting the 
surface shear forces, equation ( 10 ) reduces to 

(2+H) 

0 pU = constant ( 12 ) 


This same equation has been employed (reference 9) for estimating 
the increase in boundary— layer momentum thickness through the 
trailing— edge shock wave of an airfoil at supersonic speeds. 

It Is necessary to calculate the chordwise variation of both 
H and the momentum thickness 0 before an estimation of the growth 
of the displacement thickness 8 * may be made. As has been 
previously discussed, an analysis of the turbulent transonic 
boundary layers of reference 4 (fig. 6 ) indicated the applicability 
of the par ame ter tj for specifying the shape of the boundary layer. 
Therefore, the data of reference 4 were analyzed in terms of equa- 
tion ( 7 ) of this report. In performing this analysis, equation ( 7 ) 
was written in the form 
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.9 i(nuf)_ = 0>005 _ 0 .009 Tfj (i 3 ) 

Ur dx 

and used with the experimental values of 9 and M* to calculate the 
variation of tj. The calculated and experimental values of i) are 
compared in figure 7 . With the ohordwise variation of tj known, 
values of Hj 4 _q may be determined from equation ( 9 ) or figure 3 • 

Equation (ll) gives the correction to H^q f° r the effect of 
compressibility. Since H equals 8 */©, 5* is immediately available. 

To illustrate these results, a typical example, taken from refer- 
ence 4, is presented in figure 8 . The calculations were based on 
the experimental static pressure distribution which was converted 
to velocity distribution by using Bernoulli's equation for adiabatic 
flow. At the peak pressure, the value of H was obtained from equa— . 
tion ( 11 ) under the assumption that H^o would have been 1.4. 

This value of H was used in equation (12) to calculate the momentum- 
thickness distribution shown in figure 8 . Equation (l 3 ) then gave 
the chordwise variation of tj. Figure 3 was used to convert the 
calculated values of tj to % = o and equation ( 11 ) to determine the 

corresponding values of H. Since H equals 5*/©, 5* was easily 

determined. These values are also compared with the experimental 
measurements in figure 8 . 

The data presented by Ackeret indicate an important fact 
regarding the variation of the turbulent-boundary— layer profile 
shape through the shock wave. In the Immediate vicinity of the 
shock wave, the large pressure gradient causes the velocity profile 
to approach the shape associated with turbulent separation in low- 
speed flows. However, Immediately following this region is one 
with much smaller pressure gradients, and in this latter region the 
boundary— layer velocity profile tends to return to a shape character- 
istic of unseparated turbulent boundary layers. This tendency to 
return to a flat— plate type of velocity profile is indicated, by 
equation ( 1 3 ) and figure 7 , to be characteristic of regions of small 
pressure gradients. It thus appears possible that immediately 
behind a Buff latently intense shock wave the boundary layer may be 
separated, but re-establishment of turbulent boundary— layer flow 
will occur if the subsonic region behind the shock is sufficiently 
extensive and if the adverse pressure gradient over that region is 
small. 


DISCUSSION AND CONCLUDING REMARKS 

In the preceding analysis, a general similarity has been found 
between boundary— layer flows at low speeds and transonic speeds. 
The low— speed case considered was an airfoil section at moderate 
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angle of attack; the transonic case vas a surface over which there 
was a rather extensive region of supersonic flow. The comparison 
of the various flow regimes for these two cases will now be 
discussed. 

Laminar separation in the low— speed case is due to the 
adverse pressure gradient aft of the sharp leading— edge pressure 
peak. The potential— theory pressure distribution is so modified 
by the presence of separation that locally the pressure gradient 
is zero. The extent of the zero-gradient region is such as to 
permit the laminar flow to turn turbulent and is characterized by 
a Reynolds number run of approximately 25*000. At transonic speeds 
1 ami na-r separation is more complex because the occurrence of an 
adverse gradient and 3 ami na-r separation are two aspects of the 
same physioal phenomena; that is, it is not apparent that one is 
cause and the other effect. However, the magnitude of the flow 
deceleration immediately ahead of separation is the same as for 
low— speed flows, and in the latter case it is known that this flow 
deceleration causes IwttiI nw.-r separation. 

The examples of transonic flow considered in the present 
report have a more extensive length of laminar flow ahead of the 
transition point and hence a thicker boundary layer than the low- 
speed examples. Consequently, the transonic disturbance which 
causes 1 ami nar separation at one Reynolds number can cause abrupt 
transition at a somewhat larger Reynolds number. For moderate 
Reynolds numbers (about one million) and Mach numbers of about 1.2 
the length of the laminar— separation region is characterized by a 
Reynolds number run of about 100,000. 

At both low and transonic speeds, the boundary— layer momentum 
thickness remains constant over the constant— pressure region, and 
the displacement thickness grows in such a manner as to effectively 
change the local surface shape to one over which the pressure 
gradient would be zero. When transition occurs, there is a large 
pressure gradient and considerable pressure recovery before the 
separated flow reattaches to the surface as a turbulent boundary 
layer. Xt does not seem to have been recognized previously that 
such a large pressure recovery occurs over the transition region 
before the onset of fully developed turbulent boundary— layer flow. 
In the transition region the boundary— layer displacement thickness 
is relatively constant, while the momentum thickness increases 
markedly. 

The velocity distributions through turbulent boundary layers 
at low and transonic speeds appear to belong to the same one- 
parameter family. Thus, specifying the ratio of the local stream 
velocity to the velocity at a distance from the surface equal to 
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the local momentum thickness immediately specifies the entire 
velocity profile. This property of turbulent boundary layers 
suggests the use of the Gruschwitz method for calculating the 
chordwise variation of the velocity— prof ile shape. The available 
data indicate that such calculations give results in fair agreement 
with experimental measurements at both low and transonic speeds. 


A rough approximation determined from the Gruschwitz equation 
for the magnitude of the flow deceleration which is sufficient to 
cause an initially fully developed turbulent boundary layer to 
separate is 




senaration j _ 

Uinitial J 


It thus follows that if the boundary layer is turbulent ahead of a 
shock wave, and through the shock the local Mach number goes subsonic 
from a moderate supersonic value, then separation may exist behind 
the shock. However, the turbulent boundary layer re-establishes 
itself as unseparated flow if the adverse pressure gradient behind 
the shock is sufficiently small. This important effect of the 
pressure gradient existing behind the shock wave suggests that 
airfoil sections having small pressure gradients over the rear 
portion of the. chord will have better flow characteristics at 
transonic speeds than airfoils having large gradients in this region. 


The question arises as to whether the magnitude of the peak 
local Mach number occurring on airfoils is limited by flow separa- 
tion, The available data give no answer to this problem. However, 
the analysis indicates that, if turbulent boundary-^ layer flow 
exists ahead of the shock, the maximum deceleration through the 
shock compatible with the stability of the turbulent boundary layer 
corresponds to a decrease in Mach number from 1.5 to 1.0. Thus, if 
the flow behind the shock is to be subsonio, the properties of the 
turbulent boundary layer appear to impose a limiting value of 1.5 
for the peak Mach dumber. 

The rather rudimentary analysis of the present report cannot 
he considered as the answer to the various problems dealt with. It 
is rather intended to present a simple and coherent picture of the 
boundary— layer variation for certain airfoil configurations. The 
basic similarities between the low— speed and transonic— speed cases 
considered indicate that important simplifications and extensions 
of current boundary— layer theory are possible. 


Ames Aeronaut ioal Laboratory, 

national Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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Figure /.- Experimental pressure distribution over 
upper surface of an NAG A 66,2-216 a =0.6 
airfoil section. <2,10.1* (Data from reference 6) 
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Figure 2 - Variation of pressure coefficient and boundary -layer 
momentum and displacement thicknesses over upper surface 
of airfoil leading edge. Airfoil section , NACA 66,2-216 \ a =0.6', 
R, 900,000; Ct, iO.i T (Data from reference 6) 
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Boundary- layer shape parameter , 7/ 
Figure 3 — Variation of H with 7 . 
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Figure 4.- Experimental and calculated values of H and ff/c for a NAG A 66, 2-216, a =0.6 
oirfoif. R, 2,600,000 ; a, /O.F (Data from reference 2) 
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Figure 5. - Distribution of boundary-fayer momentum 
and displacement thicknesses and the ratio of static 
pressure to stagnation pressure through a lambda 
type shock wave. Mach number before shock wave , 
i.225 . (Data from reference 4) 
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Figure 6.— Comparison of low speed and transonic 
boundary-layer profiles. 
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(b) R e = 1478 at M= 1247. 



(c) R e = 2315 at 1.235. 


Figure 7— Var iation of boundary -layer shape parameter, 7), 
through compression shock waves. Turbulent boundary layer 
ahead of shock. (Data from reference 4) 


Mach number, 





Boundary- layer momentum and displacement thicknesses , $ and S* mm. 
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Figure 8.- Distribution of boundary-layer momentum 
and displacement thicknesses and the ratio of static 
pressure to stagnation pressure through a shock 
wave for the case of a turbulent boundary layer. 
Mach number before shock , 1279. (Data from refer- 
ence 4) 
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